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Abstract 

We employ a recent approach to the non-relativistic reduction of the electro- 
magnetic current operator in calculations of electronuclear reactions. In con- 
trast to the traditional scheme, where approximations are made for the trans- 
ferred momentum, transferred energy and initial momentum of the struck 

Q . 

5 ■ nucleon in obtaining an on-shell inspired form for the current, we treat the 

problem exactly for the transferred energy and transferred momentum. We 
calculate response functions for the reaction '^H{e, e'p)n at CEBAF (TJNAF) 
energies and find large relativistic corrections. We also show that in Plane 
Wave Impulse Approximation, it is always possible to use the full operator, 
and we present a comparison of such a limiting case with the results incor- 
porating relativistic effects to the first order in the initial momentum of the 

struck nucleon. 
25.30.Fj, 25.10+s 



X 



1 



Typeset using REVT^ 



I. INTRODUCTION 



At present, there exists a broad experimental program of electron scattering studies at 
TJNAF, MAMI, Bates, and NIKHEF aimed at understanding the short range structure of 
nuclei and the properties of nucleons in the nuclear medium. In the theoretical calculation of 
these processes, one typically performs a non-relativistic reduction of the relativistic electro- 
magnetic current operator, as the nuclear wave function is usually given in a non-relativistic 
framework. Traditionally, certain assumptions have been made for the non-relativistic re- 
duction: the momentum q transferred from the electron to the nucleus is regarded as being 
smaller than the nucleon mass mjv, and the transferred energy uo and the initial momentum 
p of the nucleon are in turn assumed to be smaller than the transferred momentum, namely, 

\q\ < tun 

U ^ < TUN (1) 

2m N 
IpI < rrtAT . 

These approximations often cannot be justified for present day experiments, since the trans- 
ferred energies and momenta may be in the GeV region and thus are comparable with or even 
bigger than the nucleon mass. This applies especially to TJNAF with its 4 GeV electron 
beam. 

In the past, there have been several attempts to improve on the reduction by expanding 
in powers oi q/niN but for a situation where also the transferred energies and not only 
the transferred momenta attain values comparable to the nucleon mass, these approaches 
are insufficient. Such a situation occurs and is indeed a common one, for instance, in near 
quasielastic kinematics. There also exist other expansion schemes 0. 

In this paper, we focus on the problem of an improved treatment of the non-relativistic 
reduction of the electromagnetic current operator. We do not discuss here, however, the 
general problem of off-shell prescriptions for the current. Instead, since we treat the non- 
relativistic reduction in a nuclear framework, we adopt the popular ansatz of using the on- 



shell form for the current and employ this form with bound and scattering wave functions 
for the nucleons in the nucleus. In so-doing this, we are, of course, using yet another off-shell 
prescription, since only with plane waves will the current matrix elements be taken strictly 
on-shell. The goal of the paper is to explore specific classes of non-relativistic approximations 
within the context of this particular ansatz, identified in this work as the "on-shell form" of 
the off-shell current. 

Some aspects of the formalism presented in this paper have been used in other recent 
publications although several approximations invoked in those studies have not been 

checked explicitly before using a realistic nuclear model. In this paper, we point out how 
to check the quality of these approximations and we present a comparison between the 
exact results (i.e. exact within the context of the model; see below for details) and the 
approximated results employing a realistic nuclear model. The calculations shown in this 
paper have been carried out for the reactions ^-ff (e, e'p)n and ^H{e, e'p)n at CEBAF energies, 
which are interesting in their own right. 

This paper is organized as follows: First, we discuss the new approach for the single- 
nucleon current that treats the problem exactly for the transferred energy and transferred 
momentum. The approach employs the on-shell form for use in off-shell calculations. In 
particular, the current is expanded in powers of the initial nucleon momentum. We show that 
for the special case of Plane Wave Impulse Approximation and exclusive (e, e'N) reactions, 
we can always employ the full form of the current operator. This allows us to check the 
quality of the expansion in powers of the initial nucleon momentum. Then we compare 
results which have been obtained using the non-relativistic reduction and the improved 
current operator, and discuss the role of the new contributions arising from the relativistic 
corrections. 
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II. THE CURRENT OPERATOR 



We start our discussion with the single-nucleon on-shell electromagnetic current operator 
and its non-relativistic reduction. Afterwards, we discuss how the current is apphed in the 
nucleus. Throughout this paper, we use the conventions of Bjorken and Drell and base 
our notation on the abbreviated treatment presented in [Q. 

It is useful to rewrite the single-nucleon current, 



J^(PA;P'A') = m(P'A') 



IrriN 



uiPA), (2) 



in a form that is more suitable for application to nuclear problems. The 4-momentum of the 
incident nucleon is = {E,p), the 4-momentum of the outgoing nucleon is P''^ = {E',p'), 
and the transferred 4-momentum is = P'^ — P^. The spin projections for incoming and 
outgoing nucleons are labeled A and A' , respectively. 

It is convenient to introduce dimensionless variables at this point: 

A = oj /2mis[ 
K, = q/2m]s[ 

r = K^-\^ (3) 
ff = p/rriN 

e = E/ttin = \Jl + 7f . 

For the outgoing nucleon, f]' and e' are defined correspondingly. Of course, these latter 
quantities can be eliminated through 4-momentum conservation. 

The Dirac and Pauli form factors are functions only of the 4-momentum transfer. Pi 2 = 
-^1,2 (t), and in the following we use the Sachs form factors Ge{t) = Ei{t) — tF2{t) and 
GM(r) = Pi(r) + P2(r). 

Also, it is useful to introduce the angle 9 between k, and fj, which allows us to obtain the 
following relations between the kinematic variables: 
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KT] COS 6 = Xe ~ T 



rie + X) 



(4) 



where 6 is defined a.s 6 = rj sin 6. For later use we also define: 

1 



1^2 



(5 + A) ^^T^ 



l+r 



(5) 
(6) 



^ — /^i 

It is our aim to obtain expressions for the single-nucleon electromagnetic current opera- 
tors J^{P; P') that occur inside the two component spin-i spinors, viz. 



(7) 



The bar over the current distinguishes an operator from its spin matrix elements. Writing 
these in the following way with an overall factor /q removed for convenience (note that J'^ in 
Eq. is a four- vector, whereas in Eq.(P) is not, although also labeled with the Lorentz 
index fi), 



(8) 

the electromagnetic current operator may then be rewritten in terms of the kinematical 
variables introduced previously: 



= (X/k) r° 



K ■ 1] 



K 



(10) 

(11) 
(12) 



+-^2 ■a){KXf])+ [{k XT])- a] 



K ■ T] 



where the coefficients (no spin dependence) and ^'^ (spin dependence) are given by: 



5 



K 



Ge + 



1 



P'iGm — -IJ'2Ge 

fJlGE + -^fl'2TGM 



I-J^^sAGm 



fi: V 2(1 + r) 
2^3 



■/il/i2G 



2K(l + r 



-/ii/i2 [Gb - G 



Ml 



(13) 



These are exact expressions for the on-shell electromagnetic current operator, also given in 
Q. So far, Eq. (|^) has only been rewritten. The conservation of the current is obvious from 
Eq. dn]). 

Later on, we will refer to the operator associated with C,o as zeroth-order charge operator, 
where the zeroth-order indicates that the operator, in this case the identity, is of zeroth order 
in 7]. Correspondingly, we call the term containing the C,o first-order spin-orbit operator, 
the term containing first-order convection current, the term containing ^[ zeroth-order 
magnetization current, the term containing ^2 first-order convective spin-orbit term, and the 
term containing ^3 second-order convective spin-orbit term. Note that all coefficients C,i and 
contain terms that are either of zeroth or first order in rj. 

Compared with the fully non-relativistic reduction discussed below (compare Eqs. (p7D ), 
there are several new types of operators in Eqs. ( p!0| - |T^ ), and the operators which are 
also present for the non-relativistic reduction, namely the zeroth-order charge operator, the 
zeroth-order magnetization current and the first-order convection current, are multiplied 
with new factors. For the charge, we now have an additional contribution, which is tradi- 
tionally called the "spin-orbit" part of the charge. It is also present in the q/rriN expansion 
schemes. In the transverse part of the current, there are two new operators with different 
spin structures (the first-order and second-order convective spin-orbit terms). 

Eqs. ( pUj ) and (|T^, which are coordinate free, can also be rewritten using a coordinate 
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system with unit vectors = k/ k, U2 = {k x if) / kt] sin 9 and -u i = ^2 x ^3 to obtain 



= 1^0 + ip'^a'^ 



uiUi + i 



"^3 



V2<^ + T^2^ ] ^2 



(14) 
(15) 



with 



K 



Ge + 



-0 tGm 



2(1 + t) 



1 



1^1 



(gm 
-rU- 



^-iGe + -^'2'tGm 



5 = 56 



2(1 + r) 

/il/i2 j-2 



(16) 



2(1 + r) 



5'\G 



M 



2 (J 



^ ( - ) Vtiii^25Gm = f^'^5^2 



In addition to expressing the current operator with respect to the coordinate system with 



unit vectors Ui,U2,'U3, the terms in Eqs. ([l^ ) , ([T5|) have also been reordered according to 
which kind of Pauh matrix they contain, so that it is easy to see which terms can interfere. 
One can read off immediately the orders of the terms that contribute to the cross section 
for electron nucleon scattering. Additionally, here it is easy to see which terms do not flip 
the spin (those involving either no Pauli matrix or a^, namely, Uo, vi, and z/2) and which do 
(those involving o"^'^, and hence cx^, namely, z/^, z/(, and z/g). 

It is instructive to evaluate the spin matrix elements of the operator given in Eqs. ([I4D 
and (|T5|) and from these to obtain the unpolarized single-nucleon responses (see also 0). 
For example, the purely longitudinal "L" response involves non-spin-flip and spin-flip ma- 
trix elements of V° squared and added incoherently (in the unpolarized responses, the two 
contributions cannot interfere): 

1 



r 



e + A)^G| + «:VGL} 



= ^{G| + m}, (17) 

(18) 

where Wi = tG\,j and W2 = [G^ + tG%j] . Similarly, the purely transverse parts of the 
current yield 

= Wi + 5'^W2 (19) 
/o [^2 + ] = = , (20) 

and from this one can see that the unpolarized "T" and "TT" responses (involving the 
sum and minus the difference of Eqs. ([T9| - p0|) , respectively) are 2Wi + 5'^W2 and —6'^W2, 
respectively. Finally, the unpolarized "TL" response is 

/2 [uou, + uX] = ^Vl + T + 6^6W2 . (21) 

Note that in Eqs. (O-PTI) the form factors Ge and Gm enter only in linear combinations of 



G|; and G\j. Terms of the type GeGm do not occur, as expected. 

When used in PWIA for the unpolarized cross section (including the special case of the 
relativistic Fermi gas model [P]) these constitute the on-shell forms for the single- nucleon 
electromagnetic response functions to be employed in concert with the nuclear spectral func- 
tion. They bear a strong resemblance to the popular off-shell prescriptions that are widely 
used in treatments of (e, e'A^) reactions |[10|-[12| and, in fact, at least in reasonably "safe" sit- 
uations such as nearly quasifree kinematics yield results that do not differ appreciably from 
those of the latter. Of course, our ultimate goal is to provide current operators and their 
matrix elements taken with interacting initial and final state wave functions, in which case 
Eqs. ([T7| - |2TD do not hold and, of course, the cross section does not factorize into single- nucleon 
responses multiplied by a spectral function. 

In electronuclear reactions, one usually assumes that the electron emits one photon which 
interacts with one of the nucleons in the nucleus. Nucleons in the nucleus are bound and 

8 



therefore off-shell; they do not fulfill the same energy-momentum relations as free nucleons. 
Currently, there exists no microscopic description of this off-shell behavior that can be 
applied for a wide range of kinematic conditions — there are only ad hoc prescriptions, 
which lead to vastly differing results for certain kinematics |TO|-[T^. 

In this paper, we want to concentrate on particular classes of non-relativistic approxima- 
tions for the current, and therefore restrict our attention to the popular ansatz of applying 
the electromagnetic current in its on-shell form. In subsequent work we intend to widen the 
scope to include other approaches to more general off-shell behaviour. 

When calculating matrix elements for the single-nucleon case, the initial and final states 
are plane waves, and therefore eigenfunctions of the operator -ff. In the nuclear case, the 
initial and final states are single-particle wave functions for nucleons in the nucleus. The 
coordinate space operator corresponding to ff is then —iV/ niiq. Unlike the plane waves 
in the case of the free nucleon, the single-particle wave functions are not eigenfunctions of 
r/ or V, and this leads to a technical problem, as the operator 17 appears several times in 
the denominators of the ^'s in Eqs. (|IB|) and the i^'s in Eqs. (|TB|). Therefore, we have 
to Taylor-expand the expressions for the current in powers of f/. Importantly, however, we 
never make expansions in k or A (and therefore r), as discussed below. 



Equations ([Tq) have been cast in forms well suited for this task. We present the results 
to first order in rj. From Eqs. @ and (|^) one finds that /^i = 1 + 0{7f) and /o = 1 + 0{rf), 



whereas one sees that /X2 = 1 + \ J j^rj cos 9 + Oijf) by using Eqs. d^-]^). We then have 



K 



v'' 



1 



Gm — l^^E 

Ge + -^tGm 



?7sin6' + 0{7f) 
7] sin 6 + 0{'rf' 



vrT7 

v^Gm + 0{r,^) (22) 



leading to the following results for the electromagnetic current operators to linear order in 





/—Ge + - 






p = 








J = 





Gm — '^Ge 



[k X r/) ■d + 0{rf') 



-i-Y- 



(k ■ a) {k X fj) 



Ge + 2^Gm 



(23) 
(24) 

(25) 



employing Eqs. ([T0|-[T2D and noting from Eq. (|) that k, = \Jt{1 + r) + rrjcosO + O^r]"^). 
Of course, when computing matrix elements of these operators and then forming bilinear 
combinations of the results to obtain the electromagnetic observables, terms of order t]"^ must 
be neglected if the operators themselves have been expanded only to order rj, since other 
terms will enter from considering the neglected 0{7f) contributions in Eqs. (p3|-p5D. 

As can be seen from Eqs. (p2|-p^, at linear order in rj we retain the spin-orbit part of 
the charge and one of the relativistic corrections to the transverse current that appeared in 
Eq. (|T2|), the first order convective spin-orbit term. 

The important point in our approach is that we have expanded only in rj, not in the 
transferred momentum k or the transferred energy A. The momentum of the initial nucleon 
will be relatively low in most cases — the dimensionless Fermi momentum rjF = pf/^n 
typically ranges from about 0.06 for deuterium to about 0.28 for heavy nuclei. However, 
there is a lot of interest in the investigation of the short-range properties of the nuclear wave 
functions which are reflected in the behavior at large momentum, and for those cases it is 
necessary to establish how good the approximation in rj is and where it does not work. 

In the next section, we will discuss the relativistic corrections to the current for the 
special case of a deuterium target. The deuteron's Fermi momentum of 55 MeV/c = 0.28 
fm~^, which corresponds to rjE ~ 0.06, is considerably smaller than the Fermi momenta of 
heavier nuclei. For many applications, the initial momenta involved are below the Fermi 
momentum, as this part of the wave function leads to large and therefore experimentally 
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more accessible cross sections. Note that the Fermi momentum itself is not a scale in the 
expansion in 77, so that when we consider heavier nuclei, the convergence at the respective 
higher Fermi momentum will be worse than the convergence at the deuteron's smaller Fermi 
momentum. For this reason, we have extended our analysis of electron scattering from the 
deuteron in the next section to rather high missing momenta in order to draw conclusions 
about the convergence of the expansion in 77 also for the region 77 > 7^|f"*eron jg 
accessible experimentally for heavier nuclei due to their larger Fermi momenta. 

Of course, for the high momentum components of the wave function, off-shell effects are 
expected also to be important and for a complete understanding of the problem have to be 
considered in addition to the relativistic corrections discussed here. 

We illustrate the importance of retaining the exact expressions for k, A and ^/T with a 
few numbers: for a momentum transfer q = 2mN, i.e. k = 1, the transferred energy under 
quasielastic conditions is a; f» 1.2 GeV corresponding to A 0.6 and r ^ 0.6. This leads 
to an extra factor n/y/r ~ 1.3 which appears in the current matrix element, and which 
may enter squared in the calculated observable, therefore leading to a 60% increase in that 
observable. Even for moderate momentum transfers, e.g. q — niN, and thus k — 0.5, the 
transferred energy is a; 400 MeV, i.e. A f« 0.2 and r ^ 0.2 and an extra factor K/^/T ^ 1.1, 
which still gives rise to a 20% increase. Apart from the factor k / ^ which originates from the 
product of the upper components and the inverse factor ^Jt j k which stems from the product 
of upper and lower components of the Dirac spinors, there also appear certain combinations 
of form factors in typical observables, namely, one often finds G% together with Gm)^- 
The latter term would be neglected in a scheme that retains only terms of leading order in 
y/r, but in fact, for example, for protons assuming dipole parameterizations for the form 
factors those two terms become equal for r ^ 0.13, which corresponds to \Q'^\ ~ (0.98 GeV/c 
)^. The importance of the relativistic corrections will become even clearer in the following 
section discussing the results for the reactions ^i?(e, e'p)n and ^H{e, e'p)n. 

As one can sec from the expressions given above, it is actually unnecessary to make 
expansions in k, Xor t at all. If one does so in spite of this, then at intermediate momentum 
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transfers the combination G'j^[ = \/^Gm should be regarded as being of leading order, and 
not of order l/mf^ as is often assumed. With these caveats under some circumstances Eq. 
(p5|) may be approximated by 

(26) 

As we will compare our results with the strict non-relativistic reduction, we quote the 
corresponding expressions: 

'^nonrel 

Jnonrel = M X (?] + Ge 

Note that the non-relativistic reduction contains both terms of zeroth order in r] and terms of 
first order in r], i.e. the convection current. The strict non-relativistic reduction is therefore 
not the lowest-order term of an expansion in 77. 

In our new approach, once the on-shell form has been adopted, the only approximation 
made is cutting off the expansion in powers of at a certain order. Of course, it is necessary 
to check how good this approximation is, e.g. if it is sufficient to use the expressions up 
to first order in rj, which are given explicitly above, or if it is necessary to include higher 
orders. One possibility to check the quality of the expansion is to evaluate the second-order 
expressions in rj and compare them with the first-order results, but as we do not know the 
convergence properties of the power series, this gives only limited information. It is clearly 
desirable to compare the first-order results with the results obtained with the full form of 
the operator. Due to the in general complicated form of the nuclear initial and final states, 
this is not possible for most cases. However, for coincidence (e, e'A^) reactions in Plane Wave 
Impulse Approximation (PWIA), it is actually possible to carry out a calculation with the 
full current operator. In this special case, one always has an outgoing plane wave for the 
knocked out nucleon in the final state, and by partially integrating the single-particle matrix 
elements of the current twice, the operator f/^ — it appears only squared, never in linear or 
cubic form — acts on the final state and thus on its eigenfunction. 
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J 



iGm X crl — G} 



K ■ 7] 



(27) 



III. RESULTS 



We have chosen to demonstrate the effects of the relativistic corrections to the current 
for the reaction ^H{e, e'p)n and ^H{e, e'p)n because of the following reasons: first, realistic 



wave functions are available in parameterized form for this nucleus |113| , ^4| and second, there 



is a lively experimental interest in the exploration of the properties of few-body systems. 



Specifically, there are several experiments planned at TJNAF with deuterium targets [|15 
However, the current operator discussed in the previous section can be applied to any nucleus 
and to any type of electronuclear reaction. 

An extensive discussion of coincidence reactions in general can be found in [0; here we 
only quote the basic formulae. The differential cross section is equal to 



rriN rUfPN 



de'dVLedVtjq j f. Sn^mi 



C" Mott free 



+h [vT'TZ^l + VTL'Tl^l^') ] , (28) 

where m,, m^v and mj are the masses of the target nucleus, the ejectile nucleon and the 
residual system, pjq and VL^ are the momentum and solid angle of the ejectile, e' is the energy 
of the detected electron and f2e is its solid angle. The helicity of the electron is denoted by 
h. The coefficients vk are the leptonic coefficients, and the TZk are the response functions 
which are defined by 

= \p{<^)f^? 

7^J;^^2J^[J;(g>,J_(g>,] 

n^f:^\Mq}M'-\J-{<i}M' 

where the J± are the spherical components of the current. 
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For our calculations, we have chosen the following kinematic conditions: the z-axis is 
parallel to g*, the missing momentum is defined as = q — p^, so that in PWIA, the 
missing momentum is equal to the negative initial momentum of the struck nucleon in the 
nucleus, pm = —p- We denote the angle between pm and g by 9, and the term "parallel kine- 
matics" indicates 6 = 0°, "perpendicular kinematics" indicates 6 = 90°, and "antiparallel 
kinematics" indicates 6 = 180°. Note that both this definition of the missing momentum 
and the definition with the other sign are used in the literature. If not stated otherwise, we 
assume that the experimental conditions are such that the kinetic energy of the outgoing 
nucleon and the angles of the missing momentum, 6 and the azimuthal angle 0, are fixed. 
If not mentioned otherwise, the kinetic energy of the outgoing proton is fixed to 1 GeV. For 
changing missing momentum, the transferred energy and momentum change accordingly — 
for the convenience of the reader, we have listed the values of the transferred energy and 
momentum for different kinematic conditions in Tables | and 0. In antiparallel kinematics, 
the transferred momentum decreases with increasing missing momentum until a point is 
reached where the kinematical limit for the process is reached, i.e. where the transferred 
4- momentum becomes timelike. This occurs at Pm = 2.7 fm~^. For this reason, the curves 
showing results for antiparallel kinematics are cut off at this value. The strong increase or 
decrease that can be observed for some cases at this point is just an artifact: when the 
responses are multiplied with the corresponding leptonic coefficients, the product goes to 
zero. 

The curves shown in this paper have been obtained using the realistic Bonn wave function 
T3| for the deuteron and we do not present results for other wave functions as we concentrate 



on the relativistic effects in the current, which are expected to be similar no matter which 
wave function one employs. In future work, other cases will be discussed. 
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A. Unpolarized Responses 



In Fig. 1^, we show the results for the longitudinal response TZl and in Fig. |^ the trans- 
verse response TZt calculated with the full current operator, the current operator to 0{ri), 
and the strict non-relativistic limit. In order to facilitate the comparison between the dif- 
ferent results in all regimes of the missing momentum, we have included both linear and 
logarithmic plots throughout this paper. For all the different kinematics, the approxima- 
tion of the electromagnetic current operator to 0{ri) is in good agreement with the results 
obtained with the full current. The two curves practically coincide up to missing momenta 
Pm ~ 1.3 — 1.5 fm^^, and the difference for higher missing momentum is small. This behavior 
is to be expected, as in PWIA the missing momentum coincides with the momentum of the 
initial nucleon inside the nucleus, which is precisely the quantity in which we have expanded 
the current operator. We have only expanded up to the first order of the initial momentum, 
and it is quite surprising that the agreement is so good up to such rather high momenta. 
This indicates that the expansion coefficients of the higher powers of t] are rather small. In 
contrast to this, the results obtained within the strict non-relativistic limit disagree strongly 
— except for parallel kinematics — with the results that include relativistic effects. From 
the linear plots, it is clear that there is a significant disagreement even for the low missing 
momenta. This disagreement of relativistic and non-relativistic treatment of the current 
increases up to almost one order of magnitude for the higher pm- 

Let us examine the validity of the assumptions stated in Eq. (^) that enter the strict 
non-relativistic reduction. For the three different kinematics presented here, the transferred 
energy uj increases with increasing missing momentum (see Table |), and therefore the as- 
sumption that UJ <^ rriN is never fulfilled; it actually evolves to > itln for higher pm- 
The initial momentum, i.e. pm, remains smaller than the nucleon mass, but comes quite 
close to it. Finally, the transferred momentum q behaves differently for the different kine- 
matic conditions (see Table H]). For antiparallel kinematics, it decreases; for perpendicular 
kinematics, it decreases slightly; and for parallel kinematics, it increases. In all cases, how- 
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ever, it remains above the nucleon mass. For parallel kinematics, in spite of the fact that 
both transferred energy and transferred momentum are larger than mjy, the relation uj < q 
holds, which in the form of ~ q'^/2mN enters the non-relativistic reduction. This is the 
reason for the comparatively better performance of the non-relativistic reduction in parallel 
kinematics. 

After these general considerations, let us turn our attention to the more specific effects 
of the relativistic corrections on the different responses. For the longitudinal response TZl, 
the relativistic corrections lead to an increase of the response function. The bulk of this 
effects stems from the factor k/a/t which is contained in and multiplies the zeroth-order 
charge operator. It reaches values of more than 1.5, and it enters squared in the response. In 
addition, a new type of operator, the first-order spin-orbit term, appears in the relativistic 
treatment. This operator is already known from the g/m^v expansions, but in our treatment 
of the relativistic effects it has a modified factor which multiplies it. Note that in a q/niN 
expansion scheme, the factor k/v^ would be treated as 1. 

For the transverse response TZt-, the relativistic corrections lead to a decrease of the 
response. This is due to the factor \fTjK that multiplies the whole transverse current. The 
fact that here \/t/ n appears rather than its inverse as in 7?.^ stems from the fact that the 
latter arises at leading order from upper x upper spinor components whereas here upper x 
lower components occur at leading order, bringing in an extra factor of r//t^. The zeroth- 
order magnetization current is by far the largest contribution to the transverse current, 
and it completely dominates TZt- The first-order convection current and the new types of 
operators which appear only in the relativistic treatment become important only for the 
interference response functions. 

In order to get a different perspective, we present in Fig. ^the longitudinal response TZl 
(left side) and the transverse response TZt (right side) for fixed momentum transfer |g| = 1.4 
GeV/c, corresponding to k = 0.75, and fixed ?/, which implies a fixed energy transfer in 
turn. For the definition of y, see |T^. The middle panels of Fig. |^ show the responses under 
quasi-elastic conditions, the top panels are kinematically "below" the quasi-elastic peak, the 
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lower panels are "above" the quasi-elastic peak. In this choice of kinematics, the angle 9 
between the missing momentum and the 2;-axis changes with the missing momentum, as 
does the kinetic energy of the outgoing proton. In contrast to the previously considered 
kinematic setting, the value of k/Vt does not change. Again, there is good agreement for 
the fully relativistic and first-order results. 

The transverse response is dominated by the zeroth-order magnetization current. There- 
fore, the difference between the non-relativistic result and the first-order result is given by 
the factor \/t/ n which has the value 0.90, 0.85, and 0.74 for the top, middle, and lower pan- 
els. Obviously, the relativistic effects are largest for the largest energy transfer. The exact 
relativistic results tend to be larger than the first-order results at higher missing momenta 
because, although the coefficient ^( of the zeroth-order magnetization current decreases with 
increasing missing momentum, the overall factor /o, which only appears in the full calcula- 
tion, increases. 

In contrast to the transverse response, the longitudinal response has significant contribu- 
tions from two terms: the zeroth-order charge operator and the first-order spin-orbit term. 
The two contributions do not interfere. The relativistic effects in the zeroth-order charge 
operator are similar to the effects in the zeroth-order magnetization current when one re- 
places the factor \/t/h with k/v^- The first-order spin-orbit term does not enter in the 
non-relativistic calculation, and its presence in the first-order results increases the difference 
of the non-relativistic and first-order treatment. 

In Fig. ^ , we show the interference response functions TZtt and TI-tl- The transverse- 
transverse response function is negative throughout the considered range, and is very small 
relative to 71l,t- As discussed in the previous section (see Eqs. (p!9| - pOD ) this behavior is 
expected, since the spin-averaged unpolarized response is proportional to 5^. As these re- 
sponses both vanish in parallel and antiparallel kinematics, we present them only for per- 
pendicular kinematics. The transverse-longitudinal response function falls between TZl,t 
and TZtt in magnitude, again as expected from the guidance provided by the spin-averaged 
unpolarized response discussed above (see Eq. (^TD ) where this contribution is seen to be 
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proportional to S. Here the full results and the results to first order in the initial nucleon 
momentum agree extremely well up to Pm = 4 fm~^. The non-relativistic results are consid- 
erably lower over the whole range of missing momenta, roughly by a factor 5. At first sight, 
this may look surprising, as one might expect that the factors which were responsible for the 
main difference between the relativistic and non-relativistic results for the longitudinal and 
the transverse responses, namely n/y/r and \fT j k would be multiplied yielding 1, and that 
therefore there would be no overall relativistic effects. However, the transverse-longitudinal 
response function has a different structure than T^.^ and T^t, as seen in Eq. (^91) . The 
response 'R-tl consists of two different contributions: one contains the product of the first- 
order spin-orbit term and the zeroth-order magnetization current, the other one contains 
the product of the zeroth-order charge operator and the first-order convection current (see 
also 0,^]). The former amounts to roughly two thirds of the total response, the latter to 
one third. As the spin-orbit operator appears only in the relativistic treatment, it is clear 
that the major contribution to 'R-tl is completely missed in the strict non-relativistic limit. 
Also, the second contribution increases when the relativistic effects are taken into account: 
besides the two factors mentioned above which cancel each other, the factor that multiplies 
the first-order convection current increases from Ge^o Ge^\t Gm-, which for protons gives 
an extra factor of approximately 1 -|- 1.4r, assuming dipole parameterization for the form 
factors. Note that those results hold under all kinematical conditions — they are not spe- 



cific for perpendicular kinematics. In |jT8|, a microscopic calculation was carried out for the 
^if(e, e'p) reaction at lower energies, and this calculation found large relativistic effects in 
T^TL- As our results show, those large relativistic effects can be attributed to the current; it 
is not necessary to incorporate relativistic dynamics in order to see large relativistic effects 



(see also |T9[] ). 

As noted above, the transverse-transverse response function is the smallest of all the 
responses. As can be seen in Fig. ^, for this response the relative difference between the 
first-order and full results is the largest one we have encountered so far. This is to be 
expected, as the transverse-transverse response is explicitly a quantity of second order in 



rj (see Eqs.(p!^-pO|)), so the second-order contributions to the current, i.e. the second-order 
convective spin-orbit term, will contribute significantly to this response, too. Up to Pm = 1-5 
fm~^, the first-order results are slightly larger than the non-relativistic results. For pm > 3 
fm~^, the first-order result starts to be smaller and disagrees more with the exact result, 
whereas the non-relativistic and exact relativistic calculations approach each other. The 
transverse-transverse response function is the product of two currents and because of this 
it has a unique structure, as seen in Eq. (|29|) . Due to this structure, the zeroth-order 
magnetization current by itself cannot contribute. This fact refiects that the response is 
of second order — it therefore cannot have a zeroth-order contribution. As the zeroth- 
order magnetization current is the largest of the components of the transverse current, this 
explains why TZj^t is so small. It also means that in the non-relativistic limit, only the 
first-order convection current can contribute. At first order in the expansion in rj, the first- 
order convective spin-orbit term (see Eq. (p4D), enters. This relativistic correction gives 
rise to a new, additive contribution to TZtt- It has the opposite sign compared with the 
first-order convection current contribution, and although it is small, at higher pm it is big 
enough to overcompensate the rise in the first-order convection current contribution due 
to relativistic effects. From the first order to the full current, the first-order convection 
current contribution rises more than the relativistic correction contribution, which in itself 
already amounts to an increase in TZtt- However, there is still another new component of 
the current that arises only in the exact expression, the second-order convective spin-orbit 
term (see Eq. (p!2|)). The contribution of this new operator alone is quite small, but due 
to the different spin structure that it has (the same as the first-order spin-orbit term in 
the charge operator), it can interfere with the zeroth-order magnetization current, which 
leads to a sizable new contribution of second order in rj. This new contribution is the 
main reason for the discrepancy between the full result and the first-order result: one of 
the major contributions to the transverse-transverse response function is present only in 
the exact treatment of the current. In the full relativistic calculation, the second-order 
convective spin-orbit term also interferes with the first-order convection current and the 
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first-order convective spin-orbit term, but those contributions are much smaller than that 
of the interference of the second-order convective spin-orbit term with the zeroth-order 
magnetization current. The fact that the non-relativistic results seem to approach the exact 
ones for higher missing momentum is therefore mere coincidence. 

In Fig. 1^, we show the interference response functions TZtt and TZtl for the same kine- 
matic conditions as used in Fig. ^. For the transverse-longitudinal response, the value of 
1^1 \fT is not significant because of the cancellation discussed earlier, and therefore the differ- 
ence between the non-relativistic and relativistic results is roughly the same for all different 
kinematics, in contrast to the longitudinal and transverse responses. For the transverse- 
transverse response, the difference between the non-relativistic and exact relativistic results 
also remains roughly constant under the different kinematic conditions, but the first-order 
results shift from the vicinity of the exact results towards the non-relativistic curve, and 
for the highest energy transfer, coincides more or less with the non-relativistic results. This 
shows that terms of second order in r/ play an important role here, namely the second-order 
convective spin-orbit term, which was discussed above. 

B. Polarized Responses 

After this discussion of the response functions that arise for unpolarized beam and unpo- 
larized target nucleus, let us discuss some of the polarization observables. As the number of 
polarization response functions is quite high, we will just discuss two representative exam- 
ples: the response 7^^/^=^ which is nonzero only for a polarized target and polarized electron 
beam, and the response TZ^^lr^. For the discussion of those responses, we will assume that 
the target is completely polarized in the Mj = 1 state, hence the superscript Mj = 1. The 
direction of the momentum transfer serves as polarization axis. 

In Fig. ^, the absolute value of the response TZ^i'~^ is shown for calculations with the 
full current operator, the current operator to 0{t]), and the strict non-relativistic limit. The 
behavior of the different curves is very similar to the case of the transverse response function 
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discussed previously: For antiparallel and perpendicular kinematics, the non-relativistic re- 
sult differs from the relativistic treatment considerably, especially for high p^- For parallel 
kinematics, all curves are quite similar, for the reasons discussed above. The agreement 
between the first-order calculation and the full calculation is excellent at lower missing 
momenta. The largest contribution to 7Zrp/~ comes from the zeroth-order magnetization 
current, and its behavior governs the behavior of the response. The relativistic treatment in- 
troduces the factor ^/T/ k, which reduces the current and the response. In the non-relativistic 
limit, the contribution from the first-order convection current is negligible. To first order in 
1], there is a very small contribution from the interference of the first-order convection current 
and the first-order convective spin-orbit term. In the exact calculation, there is an additional 
contribution from the interference between the zeroth-order magnetization current and the 
second-order convective spin-orbit term which appears only in the full treatment of the cur- 
rent operator. However, those are only small corrections to the dominant contribution of 
the zeroth-order magnetization current itself. 

For in-plane kinematics, the VJ^l~^ vanishes in parallel and antiparallel kinematics, so 
we show it only for perpendicular kinematics in Fig. |^. The most striking feature of Fig. |^ is 
the fact that in the non-relativistic limit, this response is negligible; the actual numbers are 
around 10~^^ fm^^ sr^^, whereas the response is quite sizable once the relativistic corrections 
are included. For the special case of perpendicular kinematics at the non-relativistic level, 
the only contribution to VJ^l~^ contains the product of the zeroth-order charge operator 
and the first-order convection current. For general kinematics (see the discussion of the next 
figure), the product of the zeroth-order charge operator and the zeroth-order magnetization 
current also enters. Once the first order in ri is considered, the products of the zeroth-order 
magnetization current and the first-order spin-orbit term of the charge operator and the 
product of the first-order convective spin-orbit term and the zeroth-order charge operator 
contribute to the response, leading to a magnitude that is comparable with the magnitude 
of T^TL- Again, the contribution that includes the zeroth-order magnetization current is 
the dominant one. The second contribution has the opposite sign and its size is about 20% 
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to 30% of the zeroth-order magnetization current /first-order spin-orbit contribution. For 
first order in rj, the response is negative for low Pm, it changes sign at Pm = 1-6 fm~^ and 
continues to be positive. The full result coincides with the first-order results up to pm ~ 2 
fm but it later changes sign again, therefore becoming quite different from the first-order 
result for the highest missing momenta considered here. 

In Fig. H, we show the response TZ^[r^ for fixed q and y. Under those non-perpendicular 
conditions, the non-relativistic result is of comparable size to that of the relativistic re- 
sult. The vanishing of the product of the zeroth-order charge operator and the zeroth-order 
magnetization current in perpendicular kinematics is due to the fact that it is proportional 
to spherical harmonics that vanish for 9 = 90°. Note that the spherical harmonics would 
be modified (and non-zero) once a final state interaction is included. Still, the difference 
between non-relativistic and relativistic treatment is considerable, especially for the higher 
values of y and therefore of the energy transfer. 

The current presented in this paper can be applied to heavy nuclei immediately, and 
the relativistic effects for heavy nuclei in a certain kinematical situation, i.e. for a given 
set of Pm, 0, and q and u, will be similar to what has been found for the deuteron. One 
should keep in mind that the Fermi momentum of heavy nuclei is larger than the rather 
small Fermi momentum of the deuteron. Thus, even if one takes the point of view that 
one is interested only in the region below the Fermi momentum, where the cross section is 
larger and consequently easier to access experimentally than at high pm, one will need to 
consider higher missing momenta for heavier nuclei where relativistic effects are expected to 
be larger. 

IV. CONCLUSIONS 

In this paper, we have investigated the relativistic effects that appear in the electromag- 
netic current operator. In order to perform a calculation of electronuclear reactions, we have 
chosen the popular off-shell prescription of using what we have called the "on-shell form of 
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the operator" . 

There are several calculations of particular electronuclear reactions for few-body targets 
in a microscopic fashion p!8|j20[] . Naturally, these approaches contain relativistic treatments 
of the nuclear dynamics and the current. However, it is difficult to extend these approaches 
to higher energies — the problems start above the pion emission threshold — and they 
appear unlikely to be applied to heavier nuclei in the near future. The current operator 
discussed in this paper is not restricted by these limitations, as it can be used regardless of 
energy and target. 

We have presented a formalism that allows one to perform the non-relativistic reduction 
of the free single-nucleon electromagnetic current operator without any approximation in the 
transferred momentum or transferred energy. Within the context of the chosen so-called on- 
shell form the only remaining approximation is to expand in powers of the initial momentum 
of the nucleon. In this paper, we have carried out a systematic investigation of the quality 
of this approximation by comparing the results obtained with the current to first order in 
7] with the full results for the special case of the (e, e'p) coincidence reaction in Plane Wave 
Impulse Approximation. This is the only case for which the exact relativistic calculation of 
the current can be applied to a realistic nuclear wave function. Although the assumption of 
PWIA at those energies is incomplete (see e.g. [^,^), the PWIA is a valuable and reliable 



testing ground for the relativistic effects in the operator. 

We have found very good agreement between the full relativistic treatment and the first- 
order results in almost all cases. In the special case of the PWIA, the quantity in which we 
have expanded, the momentum of the initial nucleon, coincides with the negative missing 
momentum. As expected, the 0{t]) and exact results may become somewhat different for 
missing momenta Pm ~ 3 fm and higher. 

Naturally, the size of the relativistic corrections depends on the specific kinematics. The 
non-relativistic results are far from the full calculation and even far from the first-order 
results, especially for non-parallel kinematics. 

Our results show that at GeV energies, it is necessary to take into account the relativistic 
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corrections to the electromagnetic current operator. We have also shown that it is satisfac- 
tory to include these corrections up to the first order in the initial momentum of the nucleon. 
The results in this paper have been obtained for the reactions ^H{e, e'p)n and ^H{e, e'p)n, 
but the conclusions drawn here can immediately be applied to other target nuclei and, using 
straightforward extensions of the ideas presented, to all other electronuclear reactions. 
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TABLES 

TABLE I. The transferred energy u in GeV and the dimensionless energy transfer A for 

increasing missing momentum p„i- Tlic kinetic energy of the outgoing proton is 1 GeV. The energy 
transfer does not depend on the angles of the missing momentum. 



p„/fm-i u/ GeV A 

1.00 0.53 

1 1.02 0.55 

2 1.08 0.58 

3 1.17 0.63 
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TABLE II. The transferred momentum q in GeV/c, the dimensionless momentum transfer k, 
the dimensionless negative 4-momentum transfer r, and the ratio k / -^/r for increasing missing 
momentum pm- The kinetic energy of the outgoing proton is 1 GeV. The momentum transfer de- 
pends on the angles of the missing momentum: the top part gives the values for parallel kinematics, 
the middle part for perpendicular kinematics, and the lower part for antiparallel kinematics. Note 
that for parallel kinematics, the highest accessible missing momentum is Pm = 2.7 fm~^. 



Pm/fT^ ^ q / GeV/c k t i^l^f^ 






1.70 


0.90 


0.53 


1.23 


1 


1.89 


1.01 


0.72 


1.19 


2 


2.09 


1.11 


0.91 


1.16 


3 


2.29 


1.22 


1.10 


1.16 





1.70 


0.90 


0.53 


1.23 


1 


1.68 


0.90 


0.51 


1.26 


2 


1.65 


0.88 


0.44 


1.33 


3 


1.59 


0.85 


0.33 


1.48 





1.70 


0.90 


0.53 


1.23 


1 


1.50 


0.80 


0.34 


1.37 


2 


1.30 


0.69 


0.15 


1.78 


3 


1.10 


0.59 


< 
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FIGURES 

FIG. 1. The longitudinal response function TZl is shown for parallel kinematics (a and b), 
perpendicular kinematics (c and d), and antiparallel kinematics (e and f). The solid line shows the 
response calculated with the full expression for the current operator, the dashed line is the result 
of the 0{rj) calculation and the dash-dotted line represents the result of the strict non-relativistic 
reduction. Note the different ranges of the missing momentum pm for the linear plots (left) and 
the logarithmic plots (right). 

FIG. 2. The transverse response function TZt is shown for parallel kinematics (a and b), per- 
pendicular kinematics (c and d), and antiparallel kinematics (e and f). The solid line shows the 
response calculated with the full expression for the current operator, the dashed line is the result 
of the 0{r]) calculation and the dash-dotted line represents the result of the strict non-relativistic 
reduction. Note the different ranges of the missing momcntiiiii Pm 

for the linear plots (a,c,e) and 

the logarithmic plots (b,d,f). 

FIG. 3. The longitudinal response function TZl (a,c,e) and the transverse response function TZt 
(b,d,f) are shown for fixed 3-momentum transfer |g| and different values of the y variable, which 
is defined as the negative minimal missing momentum (see e.g. [jl^). The corresponding fixed 
energy transfers (J are 0.61 GeV (a and b), 0.75 GeV (c and d), and 0.94 GeV (e and f). The 
kinematic conditions in the middle panels correspond to quasifree conditions. Note that the angle 
9 varies with changing missing momentum. The solid line shows the response calculated with the 
full expression for the current operator, the dashed line is the result of the 0{r]) calculation and 
the dash-dotted line represents the result of the strict non-relativistic reduction. 
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FIG. 4. The negative transverse-transverse response function TZtt (a and b) and the trans- 
verse-longitudinal response function TZtl (c and d) are shown for perpendicular kinematics. The 
solid line shows the response calculated with the full expression for the current operator, the dashed 
line is the result of the 0{ri) calculation and the dash-dotted line represents the result of the strict 
non-relativistic reduction. Note the different ranges of the missing momentum for the linear 
plots (a,c,e) and the logarithmic plots (b,d,f). 

FIG. 5. The negative transverse-transverse response function TZtt and the trans- 
verse-longitudinal response function TZtl are shown for fixed 3-momentum transfer \q\ and different 
values of the y variable. The corresponding fixed energy transfers u are 0.61 GeV (a and b), 0.75 
GeV (c and d), and 0.94 GeV (e and f). The kinematic conditions in the middle panels correspond 
to quasifree conditions. Note that the angle 9 varies with changing missing momentum. The solid 
line shows the response calculated with the full expression for the current operator, the dashed 
line is the result of the O (rj) calculation and the dash-dotted line represents the result of the strict 
non-relativistic reduction. 

FIG. 6. The absolute value of the response function Tl^/~^ is shown for parallel kinematics (a 
and b), perpendicular kinematics (c and d), and antiparallel kinematics (e and f). The response is 
negative for parallel and antiparallel kinematics, for perpendicular kinematics it starts out negative 
and switches sign at Pm = 1.6 fm The solid line shows the response calculated with the full 
expression for the current operator, the dashed line is the result of the 0{r]) calculation and the 
dash-dotted line represents the result of the strict non-relativistic reduction. The target is assumed 
to be completely polarized in the Mj = 1 state. Note the different ranges of the missing momentum 
Pm for the linear plots (a,c,e) and the logarithmic plots (b,d,f). 
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FIG. 7. The absolute value of the response function TZrplr is shown for perpendicular kinemat- 
ics. The solid line shows the response calculated with the full expression for the current operator, 
the dashed line is the result of the O (jj) calculation and the dash-dotted line represents the result of 
the strict non-relativistic reduction. Both first order and exact result are negative for low Pm and 
change sign at pm = 1.6 fm Whereas the first order result continues to be positive, the exact 
result changes sign again for p^ = 3.8 fm The target is assumed to be completely polarized in 
the Mj = 1 state. Note the different ranges of the missing momentum for the linear plot (a) 
and the logarithmic plot (b). 

FIG. 8. The absolute value of the response function TZ^^r^ is shown for fixed 3-momentum 
transfer \q\ and different values of the y variable. The corresponding fixed energy transfers lo are 
0.61 GeV (a), 0.75 GeV (b), and 0.94 GeV (c). The kinematic conditions in the middle panels 
correspond to quasifree conditions. Note that the angle 6 varies with changing missing momentum. 
The solid line shows the response calculated with the full expression for the current operator, the 
dashed line is the result of the 0{ri) calculation and the dash-dotted line represents the result of 
the strict non-relativistic reduction. In all cases, the response has a minus sign up to pjn = 1-5 
fm~^, and becomes and stays positive afterwards. 
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